Приложение 2 к отчету по теме "Анализ влияния параметров многофазного потока на механическую устойчивость трубопровода в условиях Новопортовского месторождения "

Динамический анализ влияния многофазного потока на устойчивость трубопровода (анализ элемента трубопровода)

Файл с анализом результатов массовых расчетов режима работы трубопровода в OLGA

версия 0.4 (Шабонас Артур, Хабибуллин Ринат) от 08.10.2017

Сравнение результатов расчета на динамическом многофазном симуляторе для двух вариантов элемента трубопровода длиной 1500 м. Рассматриваются модели трубопровода аналогичного трубопроводу Новопортовского месторождения по геометрическим параметров и параметрам флюидов. Данный файл является расчетным файлом и содержит программный код выполняющий расчет и построение графиков а также результаты расчета.

In [1]:
""" импортируем библиотеки необходимые для проведения расчетов"""
import tplread as tp
import seaborn as sns
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
%pylab inline
Populating the interactive namespace from numpy and matplotlib

Блок загрузки результатов расчетов

На динамическом симуляторе проведены параметрические расчеты. Результаты расчетов записываются в файл. Один файл соответсвует одному расчету - набору дебитов жидкости, газа и давления в элементе трубопровода. Читаем данные из каталога с результатам расчетов.

Загружается два вида моделей:

  • первая модель описывает участок трубопровода с одним "горбом" и длинным горизонтальным участком после него схема модели

  • вторая модель описывает участок трубопровода без "горба" - плоский схема модели

Обе модели содержат одинаковое количество контрольных точек в которых контролируются параметры потока. Обе модели запускались для различных расходов газа, жидкости и различных значений давления на конце трубопровода.

Задачи моделей можно сформулировать следующим образом:

  • оценить условия образования больших пробок при различных условиях работы трубопровода
  • оценить степень влияния пробок на механическую устойчивость трубопровода
  • оценить влияния "горба" на образование пробок, сравнив режимы работы в одной точке при наличии и отсутствии "горба"
In [2]:
# читаем сначала данные плоской модели
tsim = tp.TplParams('NPsimple/')
tsim.read_data()
tsim.calc_data()
# затем читаем данные модели с одним горбом модели
tpl = tp.TplParams('NPtest1/')
tpl.read_data()
tpl.calc_data()

#print(t.df_super)
160 файла найдено.
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . read done.
calc done.
160 файла найдено.
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . read done.
calc done.

Служебный блок. Проверка считанных данных.

Проверяются заголовки ключевых слов, диапазоны изменения параметров.

In [3]:
def check(t):
    # процедура вывод результатов проверки загрузки данных
    print(" файлов прочитано", len(t.files))
    print('ключевые слова для расчетов      ',t.key_list)
    print('полный список ключевых слов      ',t.key_list_all)
    print('список точек контроля параметров ',t.pipe_list)
    print('список значений давления в трубе ',t.p_end_list)
    print('список значений дебитов жидкости ',t.qliq_list)
    print('список значений дебитов газа     ',t.qgas_list)
    
check(tsim)     # проверим первый блок файлов и выведем результаты
 файлов прочитано 160
ключевые слова для расчетов       ['HOLEXP', 'USFEXP', 'USTEXP', 'USL', 'USG', 'LSBEXP', 'LSLEXP']
полный список ключевых слов       ['HOL' 'HOLEXP' 'ID' 'LSBEXP' 'LSLEXP' 'PT' 'QG' 'QGST' 'QLST' 'QOST' 'QT'
 'TM' 'TU' 'UG' 'UL' 'USFEXP' 'USG' 'USL' 'USLT' 'USTEXP']
список точек контроля параметров  ['Pipe-3' 'Pipe-5' 'Pipe-7' 'Pipe-11' 'Pipe-13' 'Pipe-17' 'Pipe-19'
 'Pipe-21' 'Pipe-23']
список значений давления в трубе  [  5.  10.]
список значений дебитов жидкости  [  7500.  10000.  12500.  15000.  17500.  20000.  22500.  25000.  27500.
  30000.]
список значений дебитов газа      [  3.   5.   7.   9.  11.  13.  15.  17.]
In [4]:
check(tpl) # проверим второй блок файлов и выведем результаты
 файлов прочитано 160
ключевые слова для расчетов       ['HOLEXP', 'USFEXP', 'USTEXP', 'USL', 'USG', 'LSBEXP', 'LSLEXP']
полный список ключевых слов       ['HOL' 'HOLEXP' 'ID' 'LSBEXP' 'LSLEXP' 'PT' 'QG' 'QGST' 'QLST' 'QOST' 'QT'
 'TM' 'TU' 'UG' 'UL' 'USFEXP' 'USG' 'USL' 'USLT' 'USTEXP']
список точек контроля параметров  ['Pipe-3' 'Pipe-5' 'Pipe-7' 'Pipe-11' 'Pipe-13' 'Pipe-17' 'Pipe-19'
 'Pipe-21' 'Pipe-23']
список значений давления в трубе  [  5.  10.]
список значений дебитов жидкости  [  7500.  10000.  12500.  15000.  17500.  20000.  22500.  25000.  27500.
  30000.]
список значений дебитов газа      [  3.   5.   7.   9.  11.  13.  15.  17.]

Соответствие приведенный названий точек контроля и их положения на трубе приведено на рисунке

схема модели

Карты влияния потока на устойчивость трубопровода

Для всех расчитанных параметров можно построить карты влияния потока на устойчивость трубопровода. Карта показывает значения коэффциента устойчивости от расхода газа и жидкости, а также отношения коэффициента для "горбатой" трубы по отношению к трубе плоской

In [5]:
def plot_map(pp, pend=10, val='mech', vm = 5, titl = ''):
    print('Контрольная точка ',pp, '===================================================================================')
    
    AA = plt.figure(figsize=(17, 6), dpi=70)
    ax1 = AA.add_subplot(121)
    
    #print('Карта для для плоского трубопровода')
    plt.title('Карта для для плоского трубопровода.' + titl)
    a=tsim.get_matr_ql_qg(pipe=pp,p_end=pend,val=val)     # расчет карты по данным моделирования
    sns.heatmap(a, annot=True,cmap="Reds", vmin=0, vmax=vm, linewidths=.5)     # построение визуального представления карты
    plt.yticks(rotation="horizontal")
    
    ax2 = AA.add_subplot(122)
    #print('Карта трубопровода с "горбом"') 
    plt.title('Карта трубопровода с "горбом".'+ titl)
    b=tpl.get_matr_ql_qg(pipe=pp,p_end=pend,val=val)      # расчет карты по данным моделирования
    sns.heatmap(b, annot=True,cmap="Greens", vmin=0, vmax=vm, linewidths=.5)    # построение визуального представления карты
    plt.yticks(rotation="horizontal")
    AA.show()
    show()
    
    AA = plt.figure(figsize=(17, 6), dpi=70)
    ax1 = AA.add_subplot(121)
    print('Карта отношений степени воздействия.'+ titl)
    plt.title('Карта отношений. '+ titl)
    c = b / a 
    sns.heatmap(c, annot=True,  cmap='Blues', vmin=0, vmax=3, linewidths=.5)
    show()
    
# вспомогательная функция для поиска номера трубы в списке
def ipipe(tpl,pp):
    pipe_ind = tsim.pipe_list.tolist().index(pp)
    return pipe_ind

def plot_trend_example(k,qg,ql,pp,p):
    AA = plt.figure(figsize=(17, 6), dpi=70)
    
    ax1 = AA.add_subplot(121)
    plt.title('труба плоская. '+ str(k))
    d1=tp.plot_trend(tsim,klist=k,qg_num=qg, ql_num=ql,pipe_num=pp,p_num=p)[4:]
    plt.plot(d1)
    plt.yticks(rotation="horizontal")
    
    ax2 = AA.add_subplot(122)
    plt.title('труба с "горбом". '+ str(k))
    d2=tp.plot_trend(tpl,klist=k,qg_num=qg, ql_num=ql,pipe_num=pp,p_num=p)[4:]
    plt.plot(d2)
    plt.yticks(rotation="horizontal")
    AA.show()
    show()
    
def plot_trend_super_example(qg,ql,pp,p):    
    AA = plt.figure(figsize=(17, 6), dpi=70)
    ax1 = AA.add_subplot(121)
    plt.title('труба плоская. Мех. влияние')
    d3 = tp.plot_trend_super(tsim,qg_num=qg, ql_num=ql,pipe_num=pp,p_num=p)
    plt.plot(d3)
    plt.yticks(rotation="horizontal")
    ax2 = AA.add_subplot(122)
    plt.title('труба с "горбом". Мех. влияние')
    d4 = tp.plot_trend_super(tpl,qg_num=qg, ql_num=ql,pipe_num=pp,p_num=p)
    plt.plot(d4)
    plt.yticks(rotation="horizontal")
    AA.show()
    show()

Построение карт влияния потока на трубопровод

Ниже строятся карты для всех контрольных точек трубопровода

In [6]:
for pp in tsim.pipe_list:           # строим для всех контрольных точек
    plot_map(pp, val='mech',titl=' Мех.влияние')
Контрольная точка  Pipe-3 ===================================================================================
Карта отношений степени воздействия. Мех.влияние
Контрольная точка  Pipe-5 ===================================================================================
Карта отношений степени воздействия. Мех.влияние
Контрольная точка  Pipe-7 ===================================================================================
Карта отношений степени воздействия. Мех.влияние
Контрольная точка  Pipe-11 ===================================================================================
Карта отношений степени воздействия. Мех.влияние
Контрольная точка  Pipe-13 ===================================================================================
Карта отношений степени воздействия. Мех.влияние
Контрольная точка  Pipe-17 ===================================================================================
Карта отношений степени воздействия. Мех.влияние
Контрольная точка  Pipe-19 ===================================================================================
Карта отношений степени воздействия. Мех.влияние
Контрольная точка  Pipe-21 ===================================================================================
Карта отношений степени воздействия. Мех.влияние
Контрольная точка  Pipe-23 ===================================================================================
Карта отношений степени воздействия. Мех.влияние

Анализ построенных карт показывает, что наибольшее влияния потока на трубу происходит на точках Pipe-7 - Pipe 13 схема модели При этом 500 метром горизонтальной трубы практически выравнивают матрицы влияния

Детальное рассмотрение отдельных вариантов

Рассмотрим детальные графики параметров для отдельных переменных, на основании которых проводился расчет карт влияния

Контрольная точка Pipe-7 (перед "горбом")

Рассмотрим карты влияния при давлении 10 атм

In [7]:
plot_map('Pipe-7',pend=10, titl=' Мех.влияние')
Контрольная точка  Pipe-7 ===================================================================================
Карта отношений степени воздействия. Мех.влияние

Видно, что непосредственно перед горбом уже наблюдается повышенная механическая нагрузка на компенсатор. Усилия выросли примерно в 2 раза (синяя карта)

Рассмотрим карты влияния при давлении 5 атм

In [8]:
plot_map('Pipe-7',pend=5, titl=' Мех.влияние')
Контрольная точка  Pipe-7 ===================================================================================
Карта отношений степени воздействия. Мех.влияние
In [ ]:
 

Картина качественно не поменялась.

Рассмотрим карты скоростей движения пробок

In [9]:
plot_map('Pipe-7',pend=10, val = 'slug_velocity', vm = 40, titl=' Скорости пробок')
Контрольная точка  Pipe-7 ===================================================================================
Карта отношений степени воздействия. Скорости пробок

Видно, что для модели скорости выше на 20% - 40%. Чем больше скорость газа - тем больше скорости движения пробок жидкости. Перед разрушением (переходом на режим расслоенного потока) скорость достигают 30 м/сек.

Рассмотрим карты "веса пробок" - разности веса пробок и пузырей между пробками

In [10]:
plot_map('Pipe-7',pend=10, val = 'slug_delta_holdup', vm = 0.3, titl=' Вес пробок')
Контрольная точка  Pipe-7 ===================================================================================
Карта отношений степени воздействия. Вес пробок

Видно, что вес пробок снижается при росте скорости газа (чем больше газа тем меньше жидкости и в пробках и в пузырях между пробками)

Построим карты средних длин пробок

In [11]:
plot_map('Pipe-7',pend=10, val = 'slug_length', vm = 200 , titl=' Длина пробок')
Контрольная точка  Pipe-7 ===================================================================================
Карта отношений степени воздействия. Длина пробок
In [12]:
plot_map('Pipe-7',pend=10, val = 'bubble_length', vm = 1000, titl = ' Расстояния между пробками')
Контрольная точка  Pipe-7 ===================================================================================
Карта отношений степени воздействия. Расстояния между пробками

Пример графиков изменения параметров потока во времени

Расчеты проводились для разных параметров потока. Виртуальное время симуляции за которое оценивались параметры потока составило 1 час. Такой подход к моделированию позволяет построить графики изменения величин во времени.

In [13]:
"""выбираем значения параметров для анализа"""
def plot_trends_t(pp_num = 0):
    qg=1           # номер значения дебита газа
    ql=8           # номер значения дебита жидкости
    pp=pp_num           # номер значения дебита участка трубы
    p=1            # номер значения дебита давления
    k=['HOL','HOLEXP']   # номер значения ключевого слова
    print('===== истинное содержание жидкости в потоке ==========')
    plot_trend_example(k,qg,ql,pp,p)
    print('===== коэффициент механического влияния на компенсатор ==========')
    plot_trend_super_example(qg,ql,pp,p)
    k=['LSLEXP','LSBEXP']   # номер значения ключевого слова
    print('===== длина пробок и расстояние между пробками ==========')
    plot_trend_example(k,qg,ql,pp,p)
    k=['USG','USL']   # номер значения ключевого слова
    print('===== скорости движения жидкости и газа ==========')
    plot_trend_example(k,qg,ql,pp,p)

    
plot_trends_t()
===== истинное содержание жидкости в потоке ==========
['HOL', 'HOLEXP'] ['Pipe-3']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['HOL', 'HOLEXP'] ['Pipe-3']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== коэффициент механического влияния на компенсатор ==========
['MECH:Pipe-3'] ['Pipe-3']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['MECH:Pipe-3'] ['Pipe-3']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== длина пробок и расстояние между пробками ==========
['LSLEXP', 'LSBEXP'] ['Pipe-3']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['LSLEXP', 'LSBEXP'] ['Pipe-3']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== скорости движения жидкости и газа ==========
['USG', 'USL'] ['Pipe-3']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['USG', 'USL'] ['Pipe-3']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
In [14]:
i=0
for pp in tsim.pipe_list:
    plot_trends_t(i)
    i=i+1
===== истинное содержание жидкости в потоке ==========
['HOL', 'HOLEXP'] ['Pipe-3']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['HOL', 'HOLEXP'] ['Pipe-3']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== коэффициент механического влияния на компенсатор ==========
['MECH:Pipe-3'] ['Pipe-3']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['MECH:Pipe-3'] ['Pipe-3']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== длина пробок и расстояние между пробками ==========
['LSLEXP', 'LSBEXP'] ['Pipe-3']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['LSLEXP', 'LSBEXP'] ['Pipe-3']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== скорости движения жидкости и газа ==========
['USG', 'USL'] ['Pipe-3']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['USG', 'USL'] ['Pipe-3']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== истинное содержание жидкости в потоке ==========
['HOL', 'HOLEXP'] ['Pipe-5']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['HOL', 'HOLEXP'] ['Pipe-5']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== коэффициент механического влияния на компенсатор ==========
['MECH:Pipe-5'] ['Pipe-5']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['MECH:Pipe-5'] ['Pipe-5']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== длина пробок и расстояние между пробками ==========
['LSLEXP', 'LSBEXP'] ['Pipe-5']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['LSLEXP', 'LSBEXP'] ['Pipe-5']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== скорости движения жидкости и газа ==========
['USG', 'USL'] ['Pipe-5']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['USG', 'USL'] ['Pipe-5']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== истинное содержание жидкости в потоке ==========
['HOL', 'HOLEXP'] ['Pipe-7']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['HOL', 'HOLEXP'] ['Pipe-7']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== коэффициент механического влияния на компенсатор ==========
['MECH:Pipe-7'] ['Pipe-7']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['MECH:Pipe-7'] ['Pipe-7']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== длина пробок и расстояние между пробками ==========
['LSLEXP', 'LSBEXP'] ['Pipe-7']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['LSLEXP', 'LSBEXP'] ['Pipe-7']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== скорости движения жидкости и газа ==========
['USG', 'USL'] ['Pipe-7']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['USG', 'USL'] ['Pipe-7']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== истинное содержание жидкости в потоке ==========
['HOL', 'HOLEXP'] ['Pipe-11']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['HOL', 'HOLEXP'] ['Pipe-11']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== коэффициент механического влияния на компенсатор ==========
['MECH:Pipe-11'] ['Pipe-11']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['MECH:Pipe-11'] ['Pipe-11']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== длина пробок и расстояние между пробками ==========
['LSLEXP', 'LSBEXP'] ['Pipe-11']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['LSLEXP', 'LSBEXP'] ['Pipe-11']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== скорости движения жидкости и газа ==========
['USG', 'USL'] ['Pipe-11']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['USG', 'USL'] ['Pipe-11']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== истинное содержание жидкости в потоке ==========
['HOL', 'HOLEXP'] ['Pipe-13']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['HOL', 'HOLEXP'] ['Pipe-13']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== коэффициент механического влияния на компенсатор ==========
['MECH:Pipe-13'] ['Pipe-13']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['MECH:Pipe-13'] ['Pipe-13']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== длина пробок и расстояние между пробками ==========
['LSLEXP', 'LSBEXP'] ['Pipe-13']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['LSLEXP', 'LSBEXP'] ['Pipe-13']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== скорости движения жидкости и газа ==========
['USG', 'USL'] ['Pipe-13']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['USG', 'USL'] ['Pipe-13']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== истинное содержание жидкости в потоке ==========
['HOL', 'HOLEXP'] ['Pipe-17']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['HOL', 'HOLEXP'] ['Pipe-17']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== коэффициент механического влияния на компенсатор ==========
['MECH:Pipe-17'] ['Pipe-17']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['MECH:Pipe-17'] ['Pipe-17']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== длина пробок и расстояние между пробками ==========
['LSLEXP', 'LSBEXP'] ['Pipe-17']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['LSLEXP', 'LSBEXP'] ['Pipe-17']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== скорости движения жидкости и газа ==========
['USG', 'USL'] ['Pipe-17']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['USG', 'USL'] ['Pipe-17']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== истинное содержание жидкости в потоке ==========
['HOL', 'HOLEXP'] ['Pipe-19']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['HOL', 'HOLEXP'] ['Pipe-19']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== коэффициент механического влияния на компенсатор ==========
['MECH:Pipe-19'] ['Pipe-19']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['MECH:Pipe-19'] ['Pipe-19']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== длина пробок и расстояние между пробками ==========
['LSLEXP', 'LSBEXP'] ['Pipe-19']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['LSLEXP', 'LSBEXP'] ['Pipe-19']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== скорости движения жидкости и газа ==========
['USG', 'USL'] ['Pipe-19']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['USG', 'USL'] ['Pipe-19']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== истинное содержание жидкости в потоке ==========
['HOL', 'HOLEXP'] ['Pipe-21']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['HOL', 'HOLEXP'] ['Pipe-21']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== коэффициент механического влияния на компенсатор ==========
['MECH:Pipe-21'] ['Pipe-21']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['MECH:Pipe-21'] ['Pipe-21']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== длина пробок и расстояние между пробками ==========
['LSLEXP', 'LSBEXP'] ['Pipe-21']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['LSLEXP', 'LSBEXP'] ['Pipe-21']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== скорости движения жидкости и газа ==========
['USG', 'USL'] ['Pipe-21']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['USG', 'USL'] ['Pipe-21']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== истинное содержание жидкости в потоке ==========
['HOL', 'HOLEXP'] ['Pipe-23']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['HOL', 'HOLEXP'] ['Pipe-23']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== коэффициент механического влияния на компенсатор ==========
['MECH:Pipe-23'] ['Pipe-23']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['MECH:Pipe-23'] ['Pipe-23']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== длина пробок и расстояние между пробками ==========
['LSLEXP', 'LSBEXP'] ['Pipe-23']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['LSLEXP', 'LSBEXP'] ['Pipe-23']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
===== скорости движения жидкости и газа ==========
['USG', 'USL'] ['Pipe-23']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
['USG', 'USL'] ['Pipe-23']  p_end =  [10.0]  q_g =  [5.0]  q_l =  [27500.0]
In [ ]:
 
In [ ]:
 
In [ ]:
 
In [ ]: